"""
Copyright (C) 2014, 申瑞珉 (Ruimin Shen)

This program is free software: you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

import sys
import numpy
import scipy.spatial


def main():
    points = numpy.loadtxt(sys.argv[1])
    cities = len(points)
    matrix = numpy.zeros([cities, cities])
    for i in range(cities):
        for j in range(cities):
            matrix[i][j] = scipy.spatial.distance.euclidean(points[i], points[j])
    numpy.savetxt(sys.argv[2], matrix, delimiter='\t')


if __name__ == '__main__':
    main()
